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Abstract: The origin of ultra-high energy (UHE) cosmic rays is still an open question. In the present 
work, we searched the possible UHE cosmic ray sources using the MAGIC telescope for the associated 
very high energy (VHE) gamma ray emission. Due to constrained propagation distance of such cosmic 
rays, we selected nearby galaxies in vicinity of the direction of the AGASA triplet and a HiRes UHE 
cosmic ray event: NGC3610 and NGC3613 (quasar remnants); Arp299 (a system of colliding galaxies). 
No significant excess in the VHE region was found found from these objects or their surrounding region. 
At multi- 100 GeV regime, the upper limits on fluxes were given against gamma ray sources in surrounding 
region. The presented limits constrain the flux of a new hypothetical source in the region, provided the 
cosmic rays are emitted from a single point-like origin. 



Introduction 

Cosmic rays up to an ^ 3 x lO'^" eV energy have 
been observed so far, while their origin remains 
unidentified. In the standard particle acceleration 
models, possible candidates are only a few types of 
the most powerful or very large-scale objects. Due 
to the lack of data on Galactic and extragalactic 
magnetic fields, it is difficult to trace back to their 
sources. However, the AGASA group claimed that 
a part of UHE cosmic ray events cluster in their 
arrival direction distribution ([13] and references 
therein). Several groups have also pointed out a 
notable angular correlation with particular class of 
objects such as BL Lacertae objects, which include 
known VHE gamma ray emitters [6]. A similar 
correlation was reported on quasar remnants [14]. 
By the quasar remnant scenario [12], the rotational 
energy of the black hole is dominantly radiated in 
MeV-TeV band whose luminosity can be a few 
orders higher than that of cosmic rays. If quasar 
remnants within several tens of Mpc are the ori- 
gin of all observed cosmic rays above 10^° eV, 
they would overshine more than any known TeV 
gamma ray source. But in the case that hundreds 
of such sources distribute over several 100 Mpc 
space, the cosmic ray flux, at least around sev- 



eral 10^^ eV, is consistent with flux limits by VHE 
gamma ray observations. 

Observation 

To find an evidence or hint of the UHE cosmic 
ray source, especially of the clustering events, we 
searched possible sources for VHE gamma rays 
with MAGIC (Major Atmospheric Gamma ray 
Imaging Cherenkov) Telescope [1]. The detec- 
tor consists of the world's largest 17-m-diameter 
reflector viewed by photo-multiplier tube camera 
with 3.5° field of view (FOV). It is located at the 
Observatorio del Roque de los Muchachos (2200 
m above sea level). La Palma, Canary Islands. 

Table 1 summarises the UHE cosmic ray clus- 
ter (AGSAA triplet) in the Plough, Ursa Major 
where with AGASA data alone three events above 
4 X 10^^ eV had been observed in a 2.5° radius. 
Recently, the HiRes group detected a 3.7 x 10^^ 
eV near the triplet and estimated the chance prob- 
ability of observing such 'quartet' to be 0.6% [9]. 
This direction coincides in the supergalactic plane. 

Along the direction of interest, we selected nearby 
objects that are capable sources of UHE cosmic ray 
observed on the Earth. To carry out delegated ob- 
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The data analysis was performed by our standard 
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analysis chain MARS (MAGIC Analysis and Re- 
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construction Software) [16]. The image of ob- 



served showers was parameterised by the conven- 



Table 1: UHE events in the cluster. Eq is the tional Hillas technique [7]. To interpret data, a 
estimated energy by each experiment in 10^^ eV number of simulated air showers were generated 
(EeV). 66 is a nominal direction error at a 68%CL. under actual telescope configurations [11]. For re- 
jection of the background hadronic shower events, 

. we used the random forest (RF) method [2]. To 

Object R.A. Decl. my D define the parameter 'hadronness' H that repre- 

NGC3610 Ilhl8m4 +58°47'0 8.5 31 sents how hadron-Uke showers look (=0 for gamma 
NGC3613 Ilhl8m6 +58°00^6 8.9 34 ray- and =1 for hadron- Uke), the RF algorithm is 
Arp299 Ilh28m5 +58 33 9 11.8 37 trained with data and simulated gamma ray shower 

samples for compatible zenith angles (29° — 45°). 
The energy of primary gamma rays was estimated 
similarly by the RF with simulated gamma ray 
showers. The energy resolution is < 25% at ener- 
gies of interest. The incoming direction of showers 
was reconstructed by the so-called DISP method 
[4]. The typical angular resolution is 0.1°. 

First we search for the VHE gamma ray emission 
directly from the observed object. To select gamma 
ray like showers, the criterion of H cut was opti- 
mised by simulated gamma ray showers by max- 
imizing their significance against surviving back- 
ground events. 6"^ distribution is compared with 
that of OFF-sources (expected background distri- 
bution) where 6 is the angular distance between 
reconstructed shower incoming direction and the 
object position. The cut for 0"^, typically ~ 0.02 
[degree^], was similarly optimised by simulated 
showers and was applied on the data. 

Figure 1 shows an example of 9^ distribution on 
Arp299 (Eest > 200 [GeV]) in which no signifi- 
cant excess was found for gamma ray signals. 
Also for other objects, no significant excesses was 
found. Therefore the upper limits of gamma ray 
fluxes were estimated for these cases. In four es- 
timated energy bins, the acceptance of gamma ray 
showers after cuts were evaluated by the simulated 
samples independent of ones used in the RF. 

The limits on the fluxes at a 95% confidence level 
(CL) are surmnarised in Table 3. The integral flux 
limits above > 200 GeV correspond to ~ 7% for 
NGC 3610 and NGC 3613 and 5% for Arp 299 to 
the Crab flux observed by MAGIC [15]. 



Table 2: Locations of the observed objects in the 
present work, my is visible magnitude. D is dis- 
tance in Mpc (Hq ~ 70 [km h"^ Mpc"^]). 



servations, we picked out the following three ob- 
jects as summarised in Table 2. 

NGC3610 is a merger remnant elliptical galaxy 
with a black hole of ^ 0.5 x 10** solar masses Mq. 
The maximum accessible cosmic ray energy is es- 
timated to be 4.4 X 10^^ eV. This galaxy has at- 
tracted attentions because of its fine structure and 
prominently warped disk, implying that a dynami- 
cal event occurred in a few 10^ years ago [8]. 
NGC36 1 3 is a quasar remnant with a black hole of 
1.6 X IO^Mq. The maximum cosmic ray energy 
is estimated to be 5.9 x 10^^ eV. This galaxy is 
located 0.7° off the most energetic triplet event. 
Arp299 (Mrk 171, VV118; [17] and references 
therein) is a starburst galaxy and a system of col- 
liding galaxies (NGC3690+IC694). It is charac- 
terised by an extreme star formation rate. The su- 
pernova rate is also as high as ~ 0.5 per year. 
Arp299 is pointed out to be a potential source to 
explain the AGASA triplet [5]. 

Operating the MAGIC telescope, we observed 
NGC3610 for 5.1 hours, NGC3613 for 5.0 hours 
each in December 2006 and Arp299 for 9.9 hours 
in January 2007. All the observations were carried 
out by the so-called wobble mode [3] in which the 
object position was tracked alternatively 0.4° east 
or west off the centre of FOV 
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Arp-299: E^>200[GeV]] 




Significance : E^^^ > 300 [GeV] 



e= [deg=l 



Figure 1: An example of 6"^ distribution (Arp299; 
-Best > 200[GeV]) . Closed circles: ON-source 
Histogram: OFF-source data, respectively. The 
vertical line: optimised cut of 9"^ for this case 
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Table 3: Upper linnits on gamma ray fluxes from 
observed objects at a 95% CL. 



To search for any emission apart from the thse ob- 
jects, the significance of the excess events by Equa- 
tion (17) in [10] was estimated for each sky point 
using the cumulative dataset. The background dis- 
tribution was modeUed by the data observed on the 
non-source region at similar zenith angles. 

Figure 2 shows the significance map for i^cst > 
300 GeV with a convolution of the angular reso- 
lution. Crosses denote positions of objects. Stars 
represent two AGASA events. The solid curve is 
a 1° circle centred on the cicumcentre of three ob- 
jects to accomodate to a moderately flat sensitiv- 
ity. Figure 3 shows the significance distribution for 
0.1° X 0.1° bins within this circle. The distribution 
is compatible with non-source hypothesis. 



Concluding remarks 

FoUowing up with AGASA-HiRes quartet detec- 
tion, we searched the possible UHE cosmic ray 




Figure 2: Significance map (.Eest > 300 [GeV]) of 
the observed region smoothed by the angular res- 
olution. ~ 0.1°. Crosses: positions of objects. 
Stars: AGASA events. The sohd curve: 1° circle 
from the cicumcentre of observed objects. 
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Constant 41 .51 ±3.04 

0.01828 + 0.06935 
Sigma 1.178+0.1 




Figure 3: The significance distribution within sohd 
curves (1° radius) in the map of Figure 2. 



sources, nearby quasar remnants NGC3610 and 
NGC3613 and a starburst galaxy Arp299, for VHE 
emission. In ~ 200 — 500 GeV energies, the upper 
limits on the gamma ray flux from each source is 
placed against each object and were ~ 8% — 12% 
Crab flux at a 95% CL. 

Over the region observed, there are no positive sig- 
nals for the VHE emission by effectively ~ 15 
hour observation. Assuming the triplet is a UHE 
cosmic ray signal from a single source, its energy 
flux yields ~ 2 eV cm"^ s"^ by the AGASA ob- 
servation. If any of the observed objects is a re- 
sponsible cosmic ray source, the present limits cor- 
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respond to ~ 3 — 5 times of the energy flux of UHE 
cosmic ray component. 

In near future, if nearby sources exist, UHE cosmic 
ray clusters will be found even clearly by higher 
quality data provided by > 1000-km^-scale ob- 
servatories. With progress of imaging Cherenkov 
telescopes as well, it is highly expected to identify 
the sources by both ways of cosmic rays physics 
and gamma ray astronomy to approach the the 
mystery of UHE cosmic ray origin. 
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